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We explore one-dimensional (1-D) samples of ultrac ord p olar molecules with attractive dipole- 
dipole interactions and show the existence of a repulsive barrier due to a strong quadrupole interac- 
tion between molecules. This barrier can stabilize a gas Sfiultracold KRb molecules and even lead 
to long-range wells supporting bound states between molecules. The properties of these wells can 
be controlled by external electric fields, allowing the forni^ion of long polymer-like chains of KRb, 
and studies of quantum phase transitions by varying tljia-^ffective interaction between molecules. 
We discuss the generalization of those results to other systems. 
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The recent achievements in the formation and manipu- 
lation of ultracold polar molecules pQ [2] have opened the 
gate to exciting new studies in several fields of physical 
sciences. Polar molecules could find uses in quantum in- 
formation [3 and precision measurements [21 , while their 
long-range and anisotropic interactions in dense sam- 
ples could provide a fertile ground for novel quantum 
gases [5j. In addition, advances in controlling the align- 
ment and orientation of polar molecules |BJ 13 enable 
the manipulation of these inter-molecular interactions, 
building a bridge between atomic, molecular, and opti- 
cal (AMO) physics, physical chemistry, and condensed 
matter physics. Until now, stable dipolar gases were 
thought to require a repulsive dipole-dipole interaction, 
such as provided by parallel dipoles perpendicular to a 
2-D plane. However, to observe interesting new correla- 
tions and phases, such as the Luttinger liquid transition 
[5] attractive interactions are needed. In this work, a 
system with such features is proposed and investigated, 
combining available techniques to produce ultracold po- 
lar molecules with the ability to precisely control their 
spatial orientation. 

In this work, we focus our attention on KRb, which has 
been trapped in relatively large amounts 1 . We first cal- 
culate the potential energy surface (PES) V^(i?, 6*1, 6*2, 0) 
of two KRb molecules approaching each other for a wide 
range of geometries. We assume that both molecules 
are in the ro-vibrational ground state of their electronic 
X^S"*" ground state, and rigid rotors, an approximation 
that is valid for i? ~ 20 a.u. (with the bond stretching by 
less than 0.15%) or larger. Fig.[l]shows the PES for three 
particular geometries when both molecular axes are in 
the same plane ((^ = 0) : the top panel depicts V when the 
molecules are aligned {9i = 6*2 = 90°), the middle panel 
for the T-orientation {9i — 0,02 = 90°), and the bottom 
panel for collinear molecules {9i — 02 — 0). Those curves 
illustrate the difference between the stronger short-range 
region where the electronic wavefunction becomes per- 
turbed and the weaker long-range region where the bond 
length of each KRb is not affected. The short-range re- 
gion is generally deep, with wells that depend strongly 



on th^^articular geometry, ranging from a few 100 K in 
Fig. [Jfiir co-planar geometries, to the tetramer K2Rb2 
boun^by ~ 4300 K with respect to the KRb-^KRb 
thresMd [5]. 

TE§<:Rb-hKRb PES was calculated at the CCSD(T) 
level theory using MOLPRO 2009.1 [lOllII], with the 
K anSJRb core electrons replaced by the Stuttgart rela- 
tivist^ECPlSSDE 52] and ECP36SDF [T^ pseudopo- 
tentia3&>respectively. The core-core and core-valence cor- 
relatTTOrt energy was modeled using a core polarization po- 
teut ifi.i ji i 'zl. Supplemental basis functions were added to 
existing basis sets for K |14j and Rb pjj. Uncontracting 
the bgsls sets, the exponents were optimized to reproduce 
the e^erimental equilibrium bond length, i?e, and dis- 
soci^j<^n energy, |16j . In the long-range region where 
the i|iteraction is small and the wave function overlap be- 
tweeC^^Dre two molecules is negligible, the interaction can 
be spiit into electrostatic and dispersion contributions. 
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FIG. ^ KRb-|-KRb PES for coplanar geometries: aligned 
(top),'^-oriented (middle), and collinear (bottom). The inset 
sketc^^ the geometry: R joins the geometric center two KRb, 
8i an^if 2 are the angles between their molecular axes and R, 
and <p*is the angle between the molecular planes. 
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Our analysis is concentrated on the coplanar geome- 
tries of Fig.[lj which depict a seemingly surprising result. 
While the top and middle panels depict the expected be- 
havior of a repulsive and slightly attractive dipole-dipole 
interaction, respectively, the collinear geometry (bottom 
panel) reveals a barrier. The existence of this barrier can 
be traced to a strong repulsive quadrupole interaction 
(see below). We also notice that it is higher (almost 7 
K in height) than that of the aligned geometry (about 4 
K). To better understand these ab initio results, we ex- 
amine the KRb+KRb interaction at large intermolecular 
separation R via the long-range expansion 
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The functions Wn may contain electrostatic {e.g. dipole 
2?, quadrupole Q, octupole O, or higher-order moments) 
and/or dispersion and induction contributions C„_i |17j . 
The first few terms (up to n = 6) are 

W^^V^ (2ciC2 - SlS2C^) , 

32? Q 



W4 



(1 + 3ciC2 - 2S1S2C0) (ci - C2) , 



W5=VO I - si S2C0 (2 - 5c? - 5c^ ) - C1C2 (6 5c? - 5c^) 

-Q'|^(l-3c?)(l-3c2)-12ciC2SiS2C^ + ^s?s2c2^|^ 
W^^C^fl + C6a(3c? + 3c2-2) + C6,2(3c?-l)(3c2-l) 

where = cosft^, Sj = sinfl^, c^^ = cos k(j). In Table |Ij 
we list the corresponding parameters obtained by least 
squares fit of the PES up to n = 8. The fitted V, Q, 
and O are also compared to ab initio values calculated at 
the all electron CCSD level of theory with the Roos ANO 
basis set [TH]. V and Q agree to better than 1%, attesting 
to the accuracy of the PES, while O is off by one order 
of magnitude, reflecting the difficulty of fitting the small 
contribution of DO compared to that of in W5 iV <^ 
Q); O does not play a significant role for KRb. Using 
Eq.Q, one can easily understand the physical origin of 
the barriers. For parallel molecules, i.e. 9i = 02 = and 
= 0, the two leading terms in V are 



ViR, 



W3 

i?3 



W5 
i?5 



(2) 



For colhnear KRb, 9 = 0, with W3 = 2V'^ and W5 = 
—6Q.^+ 4150 ~ — 6Q^, and because of the relatively 
weak V when compared to Q, the long-range attractive 
R~^ dipole interaction is overcome by a shorter-range re- 
pulsive R~^ quadrupole interaction (the attractive con- 
tribution of VO is much weaker than that of the re- 
pulsive Q^); at shorter range still, the attractive R~^ 
and higher contributions dominate (mostly due to the 
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TABLE I. Left: fit parameters (up to R~^). Right, ab initio 
values of the equihbrium separation Re, moments T>, Q, and 
O (from the geometric center). We for the collinear orienta- 
tion, and the turning points Rsr and Rq for various molecules 
AB in = of X^E+. All values are in atomic units. 



isotropic Cq^ term) and bring V down, hence the bar- 
rier. For aligned KRb, 9 = 90°, with W3 = -V^ and 
W5 = -(9/4)Q2 _|_3X)C' ~ -(9/4)2^ and afihough both 
leading contributions are repulsive, the leading repulsive 
W5 is about 3 times smaller than for the collinear case, 
hence the smaller barrier shown in Fig. [T] (the leading 
attractive Cg^o term is the same in both cases). 

Using Eq. ([I]), we study the geometries leading to a 
long-range barrier; Fig.[2]depicts its height T4op as a func- 
tion of 9i and 02 for a few twist angles (j). For (p = 0, a 
substantial barrier exists along the diagonal 9 = 61 =92, 
for small angles {9 ^ 20° or less), and for large angles 
{9 ^ 70° or more). While the barrier remains present for 
the small angle cone (^ 20°) as (j) increases, it quickly dis- 
appears for large 9. Roughly speaking, there is a barrier 
for a cone of ~ 20° for any (j), and for larger molecular 
misalignment, the barrier vanishes. A significant barrier 
can thus be maintained by aligning the molecules within 
a small angular cone, allowing ultracold KRb samples to 
remain stable and even be evaporatively cooled in vari- 
ous trap geometries (1-D when nearly collinear, and 1-D 
or 2-D when nearly aligned). 

Polar molecules can be oriented by coupling rotational 
states along a polarizing external electric field F. This 
can be achieved by using a DC electric field; however the 
small dipole moment of KRb requires field strengths that 



FIG. 2. Kop vs. 81 and 82- The main plot corresponds to a 
twist angle (f) = 0, while the two smaller plots to (j) = 40° (top) 
and 80° (bottom). Vtop is set to zero if there is no barrier. 

CD = _ 0= 40 




80 



3 



are difficult to achieve in the laboratory. An alternative 
is to add a separate polarizing laser field [19] that directly 
couples the rotational states of the molecule. Although 
this requires a much smaller DC field, non-adiabatic ef- 
fects are prominent [5], and for the sake of simplicity we 
calculate the rotational state coupling through directly 
scaling the external field. In the rigid-rotor approxima- 
tion, we get a superposition of field-free symmetric top 
states 

J,M 

labeled by their total angular momentum J with pro- 
jection M along F. After transforming the molecule- 
fixed frame potential V{R,9i,92,(t>) to the laboratory- 
fixed frame VLab(R-, ^i, ^2) [20], the field averaged poten- 
tial is found by evaluating 

F(R) = (J'M'17'|l/Lab(R,ri,f2)|JMr!). (4) 

In Fig. [3| we illustrate the effect of F on a pair of KRb 
molecules in 1-D, with 9p defined as the angle between F 
and R. For weak fields {F < 10 kV/cm), the molecules 
remain largely in their J — rotational state, the field 
mixes only small amounts of higher J states. Classically, 
they precess "wildly" on a wide cone about F, and thus 
the KRb-|-KRb interaction samples a large range of rel- 
ative angles, averaging its attractive and repulsive com- 
ponents, with the main contribution arising mainly from 
the isotropic attractive van der Waals Cefi term. This is 
depicted by the dashed lines in Fig. [sj^a) with a field of 5 
kV/cm for the aligned (left) and coUinear (right) orien- 
tations. In both cases, the interaction becomes strongly 
attractive at short distance, with the aligned geometry 
having a weak barrier (^ 1 mK) and the coUinear case 
showing no sign of a barrier. The solid lines show the 
effect of a larger electric field of 200 kV/cm; F strongly 
mixes many (~ 7) more J's and classically the molecules 
precess on a tighter cone, sampling a more restrictive 
range of angles. The anisotropic interactions do not av- 
erage to zero, and strong barriers are present for both 
aligned and coUinear cases. Fig. ^h) shows the interac- 
tion for a range oi 9p near the aligned and coUinear ori- 
entations for F = 200 kV/cm. Since the molecules then 
behave almost like rigid rods, we recover results similar 
to those in the molecular-frame. The barrier survives for 
a cone of angle 9p oi about 20° for both orientations, 
and the same conclusions about stability of 1-D and 2-D 
samples apply. For the aligned orientation, the barrier 
appears rapidly even for low fields, while larger fields 
{F > 70 kV/cm) are necessary for the coUinear case (see 
Fig life)). In both cases, the barrier grows rapidly to 
hundreds of mK, a value much higher than the typical 
kinetic energy of the trapped ultracold molecules (< 100 

mK). 

Fig. [sjb) hints at the existence of a long-range well 
for the coUinear geometry. We analyze this well in the 




FIG. 3. (a) KRb-l-KRb interaction (1-D) for weak (5 kV/cm: 
dashed lines) and strong electric fields (200 kV/cm: solid 
lines), oriented perpendicular (left) and parallel (right) to the 
intermolecular axis. The red cylinder represents the 1-D trap, 
the arrow the orientation of the field, and the sketch above 
the precessing molecules, (b) Intermolecular interaction with 
-F — 200 kV/cm for the aligned (top) and coUinear (bottom) 
geometries as a function of 9f (angle between F and R). (c) 
Height of the barrier for the aligned and coUinear orientations 
as a function of F. 



limit of infinite electric field, when the molecules are fully 
parallel {9 = 9i = 6*2 and cj) — 0) and both molecular- 
frame and laboratory-frame (for 1-D trapped molecules) 
potentials coincide {9 = 9p). We find a long-range well 
that can sustain several bound levels due to its large ex- 
tension and the large mass of the KRb molecules. For 
9 = 0, there are 7 levels, the deepest bound by nearly 2.7 
mK with inner and outer classical turning points at 110 
a.u. and 205 a.u., respectively. As 9 increases, the bar- 
rier due to the repulsion gets smaller and the well 
deepens, and the binding energies increase accordingly 
until 9 reaches 9c — 22°, at which point the barrier dis- 
appears and no more long-range bound levels exist (see 
Fig.|4f^a)). We note that for a small deviation from 9 = 0, 
the binding energies are not significantly affected, and an 
additional level v = 7 appears for 18° < 6* < 9c (inset in 
Fig.ga)). 

The variation of bound levels with 9 affects the scat- 
tering between molecules and their effective interaction. 
Assuming 9 (or 9p) as a fixed external parameter, we 
estimate the s-wave scattering phase shift S between two 
KRb, which depends on the interaction V and the wave 
number k; S < {> 0) corresponds to an effective repul- 
sive (attractive) interaction. Here, we choose k assuming 
h^k^ ^ mkBT (m: mass of KRb; kg: Boltzmann con- 
stant) for T ~ 700 nK and illustrate the effect in 
Fig. |4|^b) for the infinite field case {i.e. 9p = 9). For 
angles smaller than ^ 14.7°, the interaction is attrac- 
tive (with (5 > 0), while it becomes repulsive [5 < 0) for 
larger angles. In an ideal 1-D trap, the repulsive barrier 
at i? ~ 100 a.u. would stabilize the sample for an at- 
tractive effective interaction by preventing the molecules 
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from reaching short distances where inelastic processes 
(e.g., chemical reactions) could take place. Larger an- 
gles, but still within the stability cone, would also give 
stable samples since the effective interaction is repulsive. 
By varying the orientation of the electric field with re- 
spect to the trap axis, the behavior of the sample could be 
controlled; an effective attractive interaction would lead 
to a dense self-trapped system, i.e. a liquid-like sample, 
while an effective repulsive interaction would give a dilute 
sample behaving like a gas. Such control could probe a 
quantum phase transition between a Luttinger liquid and 
an ultracold gas [S] . One could also create a chain of KRb 
molecules weakly bound together {e.g. by using photoas- 
sociation); these would be akin to ultracold polymer-like 
chains stabilized by an external electric field and a 1-D 
trap. We note that the effective interaction can also be 
controlled by varying the magnitude of F. In Fig. |4]jc), 
we show (5 for = as a function of F for two collision 
energies corresponding to 700 nK and 900 nK, and find 
that its sign can be changed by varying F. 

Obtaining 1-D traps is challenging; assuming a har- 
monic trap in the perpendicular direction characterized 
by the frequency w, the size of the ground state a ^ 
y/h/rnio is of the order of a few 1000 a.u. for optical 
traps. Molecules at densities of 10^^ cm~^ loaded in such 
traps would be separated by roughly d ~ 1 /xm, and for 
repulsive effective interaction, the angle tan^^ ci/d 10° 
between their axes would remain within the cone of sta- 
bility. For an attractive effective interaction, the rele- 
vant angle is tan""'^ a/i?g, where Rq is the point where 
the barrier begins for two approaching molecules (see be- 
low), which requires a ^ 0.4i?g = 50 a.u. for KRb. 
Here, the sample would not be 1-D, with inelastic pro- 
cesses possibly taking place. Non-reactive species, such 
as RbCs (see below), could be considered to prevent in- 
elastic processes, or much tighter magnetic traps could be 
employed; in which case molecules in a triplet electronic 
state with a magnetic moment /i would be required. For 



KRb 



its a^E+(w = 0) state. 



60 a.u. can be 



achieved [21], and with Rq ~ 150 a.u. [22], tan ci/Rq 
would remain within the stability cone. 

The features discussed here for KRb can be general- 
ized to other polar molecules. The PES at long-range is 
well described by ([T]), and the existence of a barrier for 
perfectly coUinear molecules depends mostly on the first 
two terms (see Eq.([2])). By setting V = and neglecting 
VO, we find Rq ~ ^/3Q^JV^, the point where the R'^ 
repulsion takes over the R~^ attraction. We can also 
define a point i?sr where the shorter range R~^ attrac- 
tion takes over the R~^ repulsion, by neglecting the other 

contributions and setting V W^s/i?^ - Wq/R^ = 0, 

which gives i?sr = —Wq/W^. If Rq is outside the region 
where bonds are strongly perturbed (~ 20 a.u.) or higher 
Wn terms contribute significantly (roughly i?sr), then the 
barrier can exist. Table |T] gives Rq and i?sr for various 
polar molecules. Because Q has roughly the same am- 
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FIG. 4. (a) Long-range well energy levels vs. 9; an additional 
level V — 7 appears at 18° (inset), (b) scattering phase shift 
5 vs. for k corresponding to 900 nK for infinite F. (c) S 
for = as a function of the field strength F for collision 
energies corresponding to 700 nK and 900 nK. 



plitude for most of them, V dictates the behavior of the 
systems. Molecules with small V {e.g., LiNa and KRb) 
have a sizable Rq, and thus the existence of a barrier is 
very likely, unlike those with a large V {e.g., LiRb, LiCs 
and NaK). We also include RbCs, for which the existence 
of a barrier is uncertain. This is interesting since RbCs 
is known to not be reactive. However, a full investigation 
is required to know if a barrier exists. 

In conclusion, we found that the interaction between 
polar molecules exhibits a strong barrier when they 
are oriented about two specific geometries: aligned and 
coUinear. We also showed that the collinear setting gives 
meta-stable samples of ultracold molecules in a tight 1- 
D trap. The long-range R~^ dipolar attractive and R~^ 
quadrupolar repulsive contributions in the collinear ge- 
ometry lead to long-range wells between polar molecules 
sustaining several bound levels. Varying the orientation 
of the molecules using an external electric field allows 
for non-trivial effects, such as changing the effective in- 
teraction from repulsive to attractive, and possibly the 
phase of the sample from gas to liquid. Finally, we also 
predict the existence of the collinear barrier for various 
bi-alkali polar molecules based on the relative strength of 
the dipole and quadrupole moments. The combination of 
available techniques to produce ultracold molecules [TJ [5] 
and the ability to precisely control their spatial orienta- 
tion [F, T provide the tools to investigate such systems. 
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